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Abstract 


The  solar  wind  plasma  analyzer  on  board  Pioneer  11  provides  first 
observations  of  low-energy  positive  ions  in  the  magnetosphere  of  Saturn. 
Measurable  intensities  of  ions  within  the  energy-per-unit  charge  (E/Q) 
range  100  eV  to  8  keV  are  present  over  the  planetocentric  radial  distance 
range  y  4  R^to  16  in  the  dayside  magnetosphere.  The  plasmas  are 

found  to  be  rigidly  corotating  with  the  planet  out  to  distances  of  at 
least  10  Rg.  At  radial  distances  beyond  10  R^,  the  bulk  flows  appear  to 
be  in  the  corotation  direction  but  with  lesser  speeds  than  those  expected 
from  rigid  corotation.  At  radial  distances  beyond  the  orbit  of  Rhea  at 
8.8  Rg,  the  dominant  ions  are  most  likely  protons  and  the  corresponding 
typical  densities  and  temperatures  are  CO  .  5  cm~3  and  JLO^  ^K,  respectively, 
with  substantial  fluctuations.  Identification  of  the  raa'ss-per-unit 
charge  (M/Q)  of  the  dominant  ionspecies  is  possible  in  certain  regions 
of  Saturn's  magnetosphere  via  the  angular  distributions  of  positive  ions. 
A  large  torus  of  oxygen  ions  is  located  inside  the  orbit  of  Rhea  and  the 
densities  are  >  10  cm~3  over  the  radial  distance  range  ~  4  Rs  to  7.5  Rs. 
Density  maxima  appear  at  the  orbits  of  Dione  and  Tethys  where  oxygen 
ion  densities  are  ~  50  cm~3.  The  dominant  oxygen  charge  states  are  o2+ 
and  0^+  in  the  radial  distance  ranges  ~  4  Rs  to  7  Rs  and  7  Rs  to  8  Rs, 
respectively.  The  observations  are  suggestive  of  a  decrease  of  ion 
energies  to  values  less  than  the  instrument  energy  threshold  of  E/Q  »  100 
eV  at  the  apparent  inward  edge  of  the  torus  at  4  Rs.  Ion  temperatures 
increase  rapidly  from  ~  2  X  10$  °K  at  4  Rg  to  5  X  106  °K  at  7.3  Rs.  >  It 
is  concluded  that  the  most  likely  source  of  these  plasmas  is  the  photo- 
dissociation  of  water  frost  on  the  surface  of  the  ring  material  with  sub¬ 
sequent  ionization  of  the  products  and  radially  outward  diffusion^  The 
sources  associated  with  the  satellites  Dione  and  Tethys  are  probably  of 
lesser  strength.  >The  presence  of  this  plasma  torus  is  expected  to  have  a 
large  Influence  on  the  dynamics  of  Saturn's  magnetosphere  since  the  pres¬ 
sure  ratio  '8  of  these  plasmas  approaches  unity  at  radial  distances  as 
close  to  the  planet  as  6.5  Rg.  On  the  basis  of  these  observational  evi¬ 
dences  it  is  anticipated  that  quasi-periodic  outward  flows  of  plasma, 
accompanied  with  a  reconfiguration  of  the  magnetosphere  beyond  6.5  R^., 
will  occur  in  the  local  night  sector  in  order  to  relieve  the  plasma  pres¬ 
sure  from  accretion  of  plasma  from  the  rings. 
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I.  Introduction 


Prior  to  the  encounter  of  Saturn  with  Pioneer  11,  Jhn  situ  measure¬ 
ments  show  the  existence  of  two  planetary  magnetospheres  with  great  res¬ 
ervoirs  of  plasmas,  those  of  earth  and  Jupiter.  The  sources  and  dynamics 
of  these  plasma  regimes  of  the  two  planets  greatly  differ  and  both 
offer  exciting,  naturally-occurring  laboratories  with  an  abundance  of 
plasma  phenomena.  The  terrestrial  plasma  regimes,  the  plasma  sheet  and 
the  ring  current,  are  fed  by  both  the  solar  wind  and  the  earth's  iono¬ 
sphere.  The  dynamics  of  these  regions  are  sensitively  tuned  to  the  state 
of  the  solar  wind  which  supplies  the  energy  for  acceleration  of  these 
plasmas.  No  substantial  effect  on  the  overall  dynamics  or  as  a  plasma 
source  is  exerted  by  the  earth's  moon.  On  the  other  hand,  a  principal 
body  of  plasma  In  the  Jovian  magnetosphere  occupies  a  large  torus  In  the 
vicinity  of  the  orbit  of  the  moon  Io.  The  source  of  these  plasmas  Is  the 
sporadic  volcanic  activity  of  this  satellite.  The  ion  composition  In¬ 
cludes  those  of  sulfur  and  oxygen  as  ionization  products  of  sulfur  diox¬ 
ide.  Plasma  flows  are  dominated  by  the  corotatlonal  electric  fields  gen¬ 
erated  by  the  short  rotational  periods  and  large  magnetic  fields  associ¬ 
ated  with  the  planet.  Although  the  Ion  densities  within  the  torus  are 
high,  this  plasma  exhibits  a  relatively  low  ratio  of  plasma  and  magnetic 
field  pressures,  &  <  10“2,  as  compared  to  values  of  approximately  unity 
within  the  terrestrial  plasma  sheet  and  outer  ring  current.  With  plasma 
measurements  during  the  traversal  of  the  Saturn  magnetosphere  with 
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Pioneer  11  a  third  magnetosphere  with  a  large  reservoir  of  plasma  is 
found.  The  dominant  ion  species  appears  to  be  oxygen  and  these  plasmas 
are  located  in  a  large  torus  encircling  Saturn  and  engulfing  the  orbits 
of  Dione  and  Tethys.  Plasma  motions  are  dominated  by  rigid  corotational 
bulk  flow.  The  ions  appear  to  be  the  ionization  products  of  water  frost 
on  the  surfaces  of  the  ring  material  and  are  transported  radially  outward 
and  energized  to  form  a  large  plasma  torus.  The  dynamics  of  this  plasma 
torus  are  likely  to  differ  substantially  from  those  of  the  Io  torus  since 
the  torus  surrounding  Saturn  is  found  to  have  remarkably  high  8  values 
approaching  unity  at  radial  distances  as  small  as  ~  6  Rs.  Thus  the  pos¬ 
sibility  exists  that  the  source  of  plasma  from  photodissociation  of  water 
and  subsequent  ionization  of  the  OH  radical  by  electron  impact  with  the 
ambient  electron  velocity  distributions  exceeds  loss  rates  from  radial 
diffusion,  recombination  and  charge  exchange,  and  provides  a  plasma  pres¬ 
sure  sufficiently  large  to  disrupt  the  magnetosphere  of  Saturn.  If  so, 
Saturn's  magnetosphere  would  be  expected  to  experience  a  periodic,  per¬ 
haps  explosive  release  of  plasmas  from  the  torus.  Preliminary  results 
from  observations  of  plasmas  in  Saturn's  magnetosphere  with  Pioneer  11 
are  given  previously  by  Wolfe  ex  _al.  [1980].  We  report  here  the  results 
of  an  extensive  analysis  of  these  plasma  measurements. 
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II.  Observations 

The  measurements  of  positive  ion  velocity  distributions  within  the 
magnetosphere  of  Saturn  as  reported  here  are  gained  with  the  high- 
resolution  electrostatic  analyzer  of  the  Ames  Research  Center  solar  wind 
plasma  instrumentation  on  board  Pioneer  11.  A  description  of  this  plasma 
instrument  is  given  by  Wolfe  et  al.  [1974],  We  summarize  briefly  several 
of  the  instrument  characteristics  directly  pertinent  to  the  present  in¬ 
vestigation.  The  high-resolution  analyzer  spans  the  energy-per-unit 
charge  range  of  100  eV  ^  E/Q  £  8  keV  with  an  energy  resolution 
A(E/Q)/(E/Q>  *■  0.07  in  64  contiguous  passbands.  This  electrostatic 
analyzer  comprises  quadrispherical  plates  and  multiple  sensors  in  order 
to  provide  measurements  of  the  three-dimensional  velocity  distributions 
of  positive  ions  within  a  conical  field-of-view  with  half-angle  50.4°  and 
centered  along  the  spacecraft  spin  axis.  Two  of  the  sensors,  all 
continuous-channel  electron  multipliers,  are  designed  with  significantly 
larger  geometry  factors  relative  to  those  of  the  main  body  of  sensors 
which  are  dedicated  to  measurements  of  solar  wind  and  magnetosheath 
plasmas.  The  larger  geometry  factors  of  these  two  sensors  are  intended 
to  accommodate  the  hotter,  quasi-isotropic  plasmas  typically  found  In 
planetary  magnetospheres.  Their  f ields-of-view  are  positioned  at  the 
edges  (+  50.4°)  of  the  fan-shaped  f ield-of-view  of  the  high- 
resolution  plasma  analyzer.  The  spacecraft  spin-axis  is  directed 
parallel  to  the  central  axis  of  this  fan-shaped  field-of-view,  hence 
providing  measurements  of  the  angular  distributions  of  positive  ions 
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within  a  50.4°  cone  as  noted  above.  All  of  the  plasma  measurements 
reported  here  utilize  the  responses  of  one  of  the  sensors  with  larger 
geometry  factor  and  the  center  of  its  field-of-view  of  width  ft®  in  cone 
angle  lies  at  the  surface  of  the  effective  conical  field-of-view  of  the 
instrument.  The  trajectory  of  Pioneer  11  through  Saturn’s  dayside  mag¬ 
netosphere  and  the  orientation  of  the  sensor's  field-of-view  are  summar¬ 
ized  in  Figure  1.  Our  present  analysis  is  limited  to  the  inbound  tra¬ 
jectory  which  is  favorable  for  viewing  in  directions  near  the  corota- 
tional  flow  vector;  the  geometry  of  the  outbound  trajectory  positions  the 
field-of-view  at  large  angles,  >  110®,  to  this  bulk  flow  vector.  The 
characteristic  thermal  speeds  of  the  positive  ions  relative  to  the  bulk 
flow  speed  are  sufficiently  small  as  to  obviate  a  substantial  study  of 
plasmas  along  the  outbound  trajectory  at  radial  distances  <  10  Rs 
(Rs ,  1  Saturn  radius  51  60,000  km).  The  inbound  trajectory  as  projected 
onto  Saturn's  equatorial  plane  is  shown  in  the  left-hand  panel  of  Fig¬ 
ure  1.  Plasma  observations  for  the  period  ~  0000  to  1500  ERT  (Earth 
Received  Time  of  the  telemetry  signal)  on  September  1,  1979  are  discussed 
in  detail  here.  Definitive  ion  measurements  are  obtained  over  the  radial 
distance  range  4  Rg  to  16  Rs.  The  magnetopause  crossing  occurs  at  2209 
ERT  on  August  31  at  17.3  Rs  [Wolfe  et  d.,  1980],  The  trajectory  lies 
close  to  the  projections  of  the  Saturn-to-sun  and  Saturn-to-earth  vec¬ 
tors.  The  positions  of  four  of  Saturn's  satellites  are  also  shown  for 
times  corresponding  to  the  spacecraft's  crossing  of  their  respective 
orbits.  As  shown  in  the  right-hand  panel  of  Figure  1,  the  spacecraft 


spin  axis  is  aligned  along  the  spacecraf t-to-earth  vector.  The  axis  of 
the  f ield-of-view  of  sensor  1  is  also  shown  along  with  the  expected  coro¬ 
tation  velocity  vectors  at  5  Rs  and  15  Rs.  The  f ield-of-view  passes 
close,  <  25°,  to  this  corotation  vector  throughout  this  segment  of  the 
inbound  trajectory.  The  responses  of  the  sensor  are  accumulated  for  one- 
half  spin  period  (1  spin  period  *  7.69  seconds)  and  subsequently  tele¬ 
metered.  Corresponding  directions  of  the  axis  of  the  sensor  field-of- 
view  are  summarized  in  the  bottom  right-hand  diagram  of  Figure  1.  The 
spin  phase  angle  $  *  0“  occurs  when  this  axis  of  the  f ield-of-view  is 
parallel  to  the  ecliptic  and  ascending  toward  the  north  ecliptic  pole. 
Responses  of  the  sensor  are  accumulated  for  two  half  spins  ir/2  <  <p  < 

3x/2  and  3x/2  _<$<2*,  0<_$<  v/2,  hereafter  referred  to  as  angular 
sectors  <t>'  =  ir/2  and  *  3x/2,  respectively.  These  measurements  of 
angular  distributions  are  used  in  the  present  analysis  for  identification 
of  ion  species,  or  more  specifically,  mass-per-unit  charge,  M/Q.  The 
angular  sector  =  */2  generally  views  toward  the  corotational  bulk 
flow  vector.  For  the  measurements  of  positive  ions  in  the  dayside  mag¬ 
netosphere  the  instrument  operating  mode  constitutes  sampling  the  half¬ 
spin  sectors  <J>  *  ■  x/2  and  ■  3x/2  with  alternate  energy  passbands 
such  that  samples  of  ion  intensities  are  available  for  32  energy  pass- 
bands  for  each  of  the  two  sectors.  The  temporal  resolution  (instrument 
cycle  time)  for  gaining  each  such  series  of  measurements  at  Saturn  is 
780  seconds. 

The  responses  of  the  plasma  analyzer  (sensor  1)  to  positive  ions 
as  functions  of  energy  (E/Q)  and  time  on  September  1,  1979  are  shown  in 
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the  color  spectrograms  of  Plate  1.  These  responses  are  plotted  as  func¬ 
tions  of  the  energy  in  units  of  eV  and  the  Universal  Time  for  reception 
of  the  telemetry  signal  at  earth  (also  ERT),  and  are  color-coded  accord¬ 
ing  to  the  color  bar  on  the  right-hand  side  of  the  spectrograms.  The 
color  bar  is  calibrated  in  units  of  the  logarithm  base  10  of  the  sensor 
responses  accumulated  during  a  half-spin  of  the  spacecraft.  The  upper 
spectrogram  corresponds  to  the  spin  sector  *  x/2  (generally  viewing 
into  the  corotation  vector) ,  and  the  bottom  spectrogram  shows  responses 
for  the  sector  <f> '  »  3n/2.  The  larger  sensor  responses  for  the  spectro¬ 
gram  for  sector  »  ir/2  are  indicative  of  the  corotation  of  plasmas  in 
Saturn's  magnetosphere.  During  the  period  0000  to  0800  ERT  the  intensi¬ 
ties  of  positive  ions  are  relatively  low,  and  correspond  to  densities 
<  1  cra~3,  as  contrasted  to  those  intensities  encountered  deeper  within 
the  magnetosphere  for  the  period  ~  0800  to  1500  ERT.  The  high  responses 
after  1500  ERT  are  attributable  to  background  rates  due  to  energetic 
charged  particles.  The  instrument  is  commanded  off  for  the  period  ~ 

1630  to  2330  ERT.  Periapsis  is  located  at  1.35  Rs  and  occurs  at  1757 
ERT.  Significant  background  responses  can  also  be  seen  in  the  spectro¬ 
grams  for  the  period  ~  0900  to  1500  ERT  as  evidenced  by  the  lack  of 
modulation  of  sensor  responses  with  energy  (plate  voltage)  at  energies 
>  3  keV.  These  background  rates  are  subtracted  from  the  sensor  responses 
at  lower  energies  before  the  velocity  distribution  functions  are 
calculated.  For  examples,  the  background  counts  per  accumulation  period 
are  12  and  122  at  0943  ERT  and  1416  ERT,  respectively.  There  are 
approximately  60  complete  measurements  of  ion  energy  spectra,  with 
responses  above  background  rates,  for  the  observations  shown  in  Plate  1. 
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Analysis  of  the  Instrument  responses  proceeds  by  converting  these 
responses  to  phase  space  densities  n(V)  for  each  of  the  two  sectors  and 
by  assuming  that  M/Q  -1,  2,  4, 8  or  16.  A  Maxwellian  distribution  with 
corotational  bulk  flow  velocity  ^0,  density  N  and  temperature  T  is  em¬ 
ployed  to  find  a  best  nonlinear  least  squares  fit  for  each  of  the  above 
species.  Details  of  this  computation  are  given  in  the  Appendix.  In  gen¬ 
eral  when  observations  are  available  for  both  angular  sectors  ■  x/2 
and  -  3n/2  only  one  of  the  above  species  provides  a  good  fit.  Our 
following  discussion  presents  detailed  results  of  these  computations 
which  are  central  to  the  determinations  of  N  and  T  and  to  the  identifica¬ 
tion  of  the  species.  For  the  observations  in  the  outer  magnetosphere  for 
the  period  ~  0000  to  0800  ERT,  i.e. ,  ~  16  Rs  to  10  Rs,  few  good  fits  to  a 
corotating  plasma  of  any  species  are  acquired.  However,  all  velocity 
distributions  which  did  provide  a  good  fit,  e.g.,  at  0431  and  0602  ERT, 
are  identified  as  M/Q  ®  1,  or  H*.  It  is  likely  that  the  outer  magnetosphere 
is  not  rigidly  corotating  at  these  larger  radial  distances.  A  further 
analysis  assumes  that  the  ions  are  H4-  and  that  the  bulk  flow  VB  is  in  the 
same  direction  as  corotational  flow  but  with  a  fraction  of  the  expected 
rigid  corotation  speed.  Again  a  nonlinear  least  squares  fit  was  employed 
for  each  value  of  Vg  to  find  the  best  fit  to  the  observed  velocity  dis¬ 
tribution.  An  example  of  these  results  is  shown  in  Figure  2  for  the  vel¬ 
ocity  distribution  acquired  at  0024  ERT  (15.7  Rs).  Statistical  uncer¬ 
tainties  in  the  observed  n(V)  are  less  than  +  20Z  for  this  and  later  such 
series  of  measurements.  The  two  angular  sectors  ■  tt/2  and  ■  3*/2 
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are  noted  along  the  abscissa.  A  bulk  speed  of  0.5  VD  with  T  ■  1.2  X  10^ 
®K  and  N  =*  0.19  cm-^  provides  a  good  fit  to  the  measurements  If  the 
high-speed  tails  of  the  distribution  at  V  >  3  X  10?  cm/sec  are  excluded. 
On  the  other  hand,  not  all  of  the  observed  velocity  distributions  In  the 
outer  magnetosphere  can  be  fit  with  a  Maxwellian  distribution  flowing 
with  a  fraction  of  the  corotation  speed.  An  example  of  a  poor  fit  is 
shown  for  0117  ERT  in  Figure  3.  Each  velocity  distribution  provides  a 
reasonable  fit  to  the  angular  sector  <t>'  =«  x/2,  excluding  the  high-speed 
tail,  but  no  apparent  fit  for  the  complementary  sector  =  3x/2  is 
found.  There  are  three  immediate  possibilities  for  the  origins  of  this 
poor  fit:  temporal  variations,  a  non-Maxwe Ilian  velocity  distribution  or 
missing  telemetry  for  the  sector  4> *  *  3ir/2  in  the  velocity  range  ~  1.8 
to  2.8  X  10?  cm/sec.  However,  approximately  40%  of  the  observed  velocity 
distributions  with  usable  measurements  for  both  angular  sectors  can  be 
fit  with  such  a  Maxwellian  with  bulk  speed  Vg  and  directed  along  the 
corotational  vector.  This  flow  speed  ranged  from  0.3  to  0.8  VQ.  Rigid 
corotation  of  a  hydrogen  plasma  is  observed  later  during  the  period  0800 
to  1010  ERT  which  corresponds  to  a  radial  distance  range  of  10  Rs  to  8.2 
Rs.  A  substantial  fraction  of  these  velocity  distributions  are  well  fit 
with  the  computed  distributions  for  Vg  -  VQ.  Thus  the  plasmas  within 
Saturn' 8  magnetosphere  are  rigidly  corotating  with  the  planet  at  radial 
distances  <  10  Rg.  Beyond  10  to  12  Rs  the  plasma  bulk  velocities  appear 
to  be  along  the  corotation  direction  but  with  a  fraction,  0.3  to  0.8,  of 
the  expected  corotation  speed. 
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For  the  period  after  the  encounter  with  corotating  hydrogen  plasmas, 
i.e.,  after  1010  ERT,  heavy  ions  and  increasing  densities  are  detected 
with  the  plasma  instrument  (see  Plate  1).  Two  zones  which  are  each 
distinguished  by  differing  M/Q  values  for  the  ions  are  found.  These  two 
zones  span  the  periods  1010  to  1130  ERT  (8.2  Rs  to  7.1  Rs)  and  1130  to  1442 
ERT  (7.1  Rs  to  4.1  Rs),  respectively.  Examples  for  the  best  fits  to  the 
observed  velocity  distributions  for  various  assumed  M/Q  values  are  shown 
in  Figures  4  and  5.  Best  fits  for  four  ion  species  for  the  ion  velocity 
distributions  observed  at  1127  ERT  (7.1  Rs)  are  shown  in  Figure  4.  It 
should  be  noted  that  the  two  velocity  distributions  corresponding  to 
sectors  =  x/2  and  =  3x/2  provide  a  very  sensitive  measure  of  the 
M/Q  of  the  ion  species.  The  density  N  and  temperature  T  for  each  best 
fit  is  also  given  in  Figure  4.  Poor  fits  to  the  observed  velocity 
distributions  are  yielded  by  M/Q  values  of  1  (H-*-)  and  16  (0+)  whereas 
considerably  better,  but  not  perfect  fits  are  provided  by  M/Q  values  of  4 
(He+)  and  8  (0^+).  Inspection  of  these  latter  two  fits  yields  an 
assessment  of  the  M/Q  of  the  dominant  species  as  5  to  6  amu/unit  charge 
with  a  density  N  =  10  cm- 3  and  temperature  T  2  1.5  X  10^  °K.  An  example 
of  similar  calculations  for  the  second  ion  zone  at  lesser  radial 
distances  is  given  in  Figure  5.  This  velocity  distribution  at  1219  ERT 
(6.4  Rs)  is  very  well  fit  with  an  ion  distribution  with  M/Q  =  8,  N  ** 

40.8  cra“3  and  x  =  1.6  x  10^  °K  with  the  exception  of  the  high-speed  tails 
at  V  y  1.2  X  10^  cm/sec.  We  note  here  the  known  presence  of  water  frost 

/w 

on  the  satellites  Dione  and  Tethys  [Morrison  ££  al_.  ,  1976;  Fink  et  al. , 
1976]  which  are  embedded  in  this  large  plasma  torus,  and  on  the  nearby 


ring  material  (Pilcher  al_. ,  1*570;  Kuiper  £t  aK  ,  19701  ,  and  identify 
the  ions  with  M/Q  a  5-6  as  0^+  (M/Q  *  5.3)  and  those  with  M/Q  =  8  as  0^+. 
For  the  time  period  1245  to  1442  ERT  (6.0  Rs  to  4.1  Rs),  the  plasma  has 
cooled  to  the  extent  that  velocity  distributions  cannot  be  determined  for 
the  angular  sector  *  3n/2  (cf.  Plate  1).  Without  this  complementary 
angular  measurement,  the  identification  of  ion  species  becomes  less 
certain  since  it  relies  upon  the  velocity  distributions  for  sector  * 
x/2  only.  An  example  of  observations  of  ion  velocity  distributions  in 
these  cooler  plasmas  is  shown  in  Figure  6,  together  with  best  fits  to 
corotating  Maxwellian  distributions  for  several  ions.  The  fit  for  0+  is 
poor  since  no  peak  in  the  velocity  distribution  is  observed  and  hence 
ions  with  M/Q  y  16  can  be  eliminated  as  possibilities.  Reasonably 
good  fits  are  obtained  for  the  M/Q  range  ~  4  to  8  but  those  for  M/Q  -  1 
and  2  are  poorer  since  a  considerable  fraction  of  the  observed  densities 
fall  below  the  best-fit  Maxwellian  distributions.  These  comparisons  are 
consistent,  but  not  exclusively  supportive,  of  our  conclusion  that  the 
ions  within  the  radial  distance  range  7.1  Rs  to  4.1  Rs  are  0^+.  At 
lesser  radial  distances  than  4.1  Rs,  no  responses  of  the  plasma  instru¬ 
ment  to  positive  ions  in  the  energy  range  100  eV  _<  E/Q  <_  8  keV  are 
discernible  (cf.  spectrograms  of  Plate  1  for  times  later  than  1442  ERT). 
However,  the  rapidly  decreasing  average  ion  energy  with  decreasing 
radial  distance  for  this  time  period  as  clearly  evidenced  in  the  top 
spectrogram  of  Plate  1  is  suggestive  of  a  decrease  of  the  ion  average 
energy  to  values  below  the  threshold  of  the  instrument  at  100  eV-per- 


unit  charge 


We  note  again  the  substantial  high-speed  tails  of  several  of  the 
velocity  distributions  presented  here,  for  example,  those  for  the  velo¬ 
city  distribution  at  V  )  1.2  X  10?  cm/sec  in  the  lower  left-hand  panel  of 
Figure  5.  The  dominant  ion  species  is  f)2+  in  this  region.  On  the 
other  hand,  the  densities  within  the  high-speed  tails  are  sufficiently 
large  to  encourage  an  attempt  at  identification  of  the  associated  ion 
species.  The  velocity  distribution  observed  at  1206  ERT  (6.6  R8)  is 
chosen  here  for  analysis.  The  dominant  ion  is  also  with  density  and 
temperature,  26.8  cm~3  and  1.6  X  10^  °K,  respectively.  The  results  of 
calculations  for  best-fit  corotating  Maxwellians  for  several  ion  species 
are  summarized  in  Figure  7.  No  reasonable  fits  are  obtained  for  M  *  1,  2 
or  16  amu  thus  eliminating  H"*"  and  0+  as  the  secondary  ion  component  that 
appears  in  the  velocity  distribution.  The  best  fits  for  M  =  4  and  8  amu 
are  both  sufficiently  good  that  0^+  (M/Q  **  8)  and  0^+  (M/Q  ■  5.3)  cannot 
be  distinguished,  although  M  ■  5-6  fit  (not  shown)  provides  a  somewhat 
better  fit  than  M  *  4  or  8.  The  corresponding  temperatures  and  densities 
are  4.8  X  10^  °K  and  6.9  cm“^,  respectively,  if  the  ions  are  0^+,  and  8.1 
X  106  ®K  and  6.9  cra~3  for  0^+.  (See  Appendix  for  treatment  of  cases  for 
Q  >  1.)  Thus  the  second  ion  component  is  substantial  in  density  relative 
to  that  of  the  0^+  for  this  example  and  is  probably  either  a  hotter  velo¬ 
city  distribution  of  0^+  or  a  coexisting  0^+  ion  plasma.  A  reasonable 
upper  bound  to  the  H+  densities  within  the  energy  range  100  eV  to  8  keV 
at  this  position  in  the  plasma  torus  is  ~  5  cm--*. 
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The  Individual  measurements  of  positive  Ion  spectra,  a  total  of  ap¬ 
proximately  60  such  observations,  were  each  subjected  to  the  same  exten¬ 
sive  nonlinear  least  square  fit  as  that  for  the  examples  given  above. 
These  examples  qualitatively  cover  the  entire  range  of  types  of  posltlve- 
lon  velocity  distributions  that  are  observed  In  the  dayside  magneto¬ 
sphere.  A  summary  of  the  complete  series  of  observations  Is  given  In 
Figure  8  which  displays  ion  density,  dominant  ion  species  and  Ion  temper¬ 
ature  as  functions  of  planetocentric  radial  distance.  Ion  densities  over 
the  radial  distances  ~  9  to  16  Rs,  and  outside  the  orbit  of  Rhea,  gen¬ 
erally  ranged  from  ~  0.2  to  0,5  cm--*.  These  densities  raonotonlcally  In¬ 
crease  from  0.6  cm”  3  at  Rhea's  orbit  to  45  cm-^  at  the  orbit  of  Dlone. 
Densities  are  ~  10  to  50  cm“^  in  the  plasma  torus  at  4.2  to  7.5  Rs.  Two 
density  maxima  appear  to  be  associated  with  Tethys  and  Dlone.  The  Ion 
density  decreases  rapidly  with  decreasing  radial  distance  to  Instrument 
threshold  values  just  outside  the  orbit  of  Enceladus.  The  center  panel 
of  Figure  8  shows  our  best  evaluation  of  ion  species  as  per  the  proce¬ 
dures  discussed  above.  Periods  for  which  there  are  (1)  angular  distribu¬ 
tions  with  good  fits  to  a  corotating  Maxwellian  distribution,  (2)  angular 
distributions  with  poor  fits  to  rigid  corotation  and  (3)  no  available 
angular  distributions  are  specifically  designated  In  Figure  8.  Our  in¬ 
terpretation  of  the  dominant  Ion  species  is  also  shown:  0^+  at  4.1  to 
7.1  R8,  o3+  at  7.1  to  8.2  R8,  and  ff*1  at  8,2  to  16  Rs,  subject  to  the 
assumptions  of  our  analysis  as  presented  above  for  specific  examples. 

Ion  temperatures  outside  the  orbit  of  Rhea  fluctuate  greatly  and  range 
from  ~  3  X  10^  to  5  X  10^  CK.  A  maximum  in  the  ion  temperature  of  ~  6  X 


106  °K  is  centered  at  9.8  Ra.  Inside  the  orbit  of  Rhea,  ion  temperatures 
monotonically  increase  to  a  maximum  temperature  of  5  X  10^  °K  at  7.3  RB 
and  subsequently  decrease  with  decreasing  radial  distance  to  ~  2.5  X  lO1* 
#K  at  4.1  Rs  just  outside  the  orbit  of  Enceladus.  The  ion  temperatures 
in  this  latter  region  at  4.1  to  7.3  Rs  exhibit  an  extremely  steep  grad¬ 
ient  and  are  roughly  proportional  to  R^.  Overall  an  inspection  of  Figure 
8  shows  that  the  major  plasma  feature  of  the  dayside  magnetosphere  is 
the  large  plasma  torus  of  heavy  ions,  dominantly  o2+  and  0^+,  at  ~  4  to  8 
Rs . 

Other  remarkable  features  of  the  plasma  torus  are  the  magnitudes  of 
the  ion  energy  densities,  e,  and  the  ratios  of  plasma  and  magnetic  field 
pressures,  8,  at  radial  distances  deep  within  Saturn's  magnetosphere. 
These  parameters  are  shown  as  functions  of  radial  distance  in  Figure  9. 
The  scale  for  e,  in  units  of  ergs/cra^,  is  given  along  the  left-hand  ordi¬ 
nate  and  for  8,  along  the  right-hand  ordinate.  Energy  densities  beyond 
the  orbit  of  Rhea  range  from  ~  2  X  10-**  to  2  X  10“^®  erg/cm^.  The  cor¬ 
responding  8-values  are  0.02  to  0.5.  The  pressure  ratio  8  has  been  cal- 

3 

culated  using  a  centered  dipole  field  with  magnetic  moment  0.2  gauss-Rs 
[cf.  Smith  et  al. ,  1980].  The  actual  fields  at  the  larger  radial  dis¬ 
tances  near  the  magnetopause  are  somewhat  greater  than  those  calculated 
with  this  dipole  field  approximation,  but  will  not  qualitatively  affect 
the  relatively  high  values  for  8.  There  is  a  relative  minimum  in  the 
values  for  8,  0.02,  at  8.8  to  9.7  Rs  just  outside  the  orbit  of  Rhea. 
Maximum  energy  densities  are  positioned  at  ~  6.0  to  7.9  Rs  within  the 
plasma  torus  and  attain  values  as  high  as  1.1  X  10“^  ergs/cm^  at  Dione's 
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orbit.  There  is  a  more  modest  maximum  of  energy  density,  ~  4  X  10“^ 
erg/cm3,  near  the  orbit  of  Tethys.  The  values  for  8  are  particularly 
high,  >  0.2,  over  the  radial  distance  range  ~  6.0  to  8.0  Rs  and  approach 
unity  in  the  vicinity  of  7.5  Rg.  To  this  stress  on  the  magnetic  field 
must  be  added  that  due  to  the  rapid  corotation  of  the  plasma.  It  appears 
likely  that,  in  view  of  these  large  values  for  3  and  the  nature  of  the 
mechanism  that  provides  a  continuous  source  of  ions  to  the  plasma  torus, 
the  magnetosphere  of  Saturn  is  marginally  stable  or  unstable  for  confine¬ 
ment  of  torus  plasmas  at  radial  distances  as  close  to  the  planet  as  6.5  Rg. 


P 


III.  Discussion 


First  measurements  of  plasmas,  specifically  low-energy  positive  ions, 
within  Saturn's  magnetosphere  are  reported.  The  observations  are  gained 
with  an  electrostatic  analyzer  that  is  intended  for  comprehensive  deter¬ 
minations  of  the  velocity  distributions  of  solar  wind  ions  and  electrons. 
Positive  ion  intensities  within  the  dayside  magnetosphere  are  sufficiently 
large  that  clear,  usable  analyzer  responses  to  these  ions  occur  over  a 
radial  distance  range  of  ~  4  to  16  Rs.  An  electrostatic  analyzer  is  in¬ 
capable  of  determining  the  ion  species:  however,  the  observed  corotation  of 
the  plasma  allows  an  assessment  of  mass-per-unit-charge,  M/0,  of  the  dom¬ 
inant  ion  species.  The  plasma  is  observed  to  rigidly  corotate  with  the 
planet  to  planetocentric  radial  distances  of  ~  10  Rs.  Beyond  10  Rg  the 
bulk  flow  speed  varies  and  appears  to  be  ~  0.3  to  0.8  of  the  expected  rigid 
corotation  speed  V0,  presumably  due  to  viscous  drag  from  the  solar  wind 
and/or  mass  loading  from  the  magnetospheric  plasmas.  The  dominant  ion 
species  beyond  a  radial  distance  of  ~  8  Rs  appears  to  be  H+,  subject  to  the 
conditions  of  the  interpretation  that  are  given  In  the  previous  discussion 
of  results.  Inside  radial  distances  of  8  Rs,  heavy  Ions  with  M/Q  ~  8  and 
5-6,  respectively,  are  detected  within  a  large,  relatively  dense  plasma 
torus  surrounding  Saturn.  This  plasma  torus  extends  from  ~  4  Rg  to  9  Rs 
and  Is  approximately  bounded  by  the  orbits  of  Enceladus  and  Rhea.  Peak  ion 
densities  within  the  plasma  torus  are  ~  50  cm-^,  whereas  beyond  the  torus 
in  the  outer  magnetosphere  these  densities  are  typically  less  by  a  factor 
of  100,  or  ~  0.5  cm-^.  A  maximum  in  ion  temperature  of  ~  5  X  10®  °K  in 
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the  torus  Is  located  at  ~  7.5  Rs  and  a  steep  gradient  of  ion  temperatures 
occurs  with  decreasing  radial  distance  to  4.1  Rs  where  the  temperature  Is 
2.5  X  10^  °K.  The  ratios  of  plasma  pressure  to  magnetic  field  pressure, 
8,  are  remarkably  high  and  are  near  values  of  unity  as  deep  within 
Saturn's  magnetosphere  as  6.5  Rs.  This  situation  appears  to  indicate 
that  Saturn's  magnetospherlc  configuration  is  marginally  stable  or  un¬ 
stable  at  these  radial  distances  and  beyond.  It  Is  suggested  that  this 
large  plasma  pressure  Is  periodically  relieved,  perhaps  with  a  dramatic 
outflow  of  plasma,  in  the  tailward  sector  of  the  magnetosphere  where 
magnetic  stresses  are  less  likely  to  contain  these  plasmas.  Two  distinct 
maxima  of  Ion  densities  occur  within  the  plasma  torus,  one  each  at  the 
orbits  of  Dione  and  Tethys.  In  lieu  of  another  known  source  of  heavy 
ions  in  this  region,  with  the  possible  exception  of  the  E  ring,  we  con¬ 
clude  that  the  source  for  the  plasma  torus  is  Saturn's  rings  or  the 
satellites  Dione  and  Tethys. 

There  are  two  principal  measurements  of  other  phenomena  In  the  day- 
side  magnetosphere  with  Pioneer  11,  that  are  directly  correlated  at  first 
Inspection  with  the  presence  of  the  plasma  torus:  a  depression  In  mag¬ 
netic  field  magnitudes  below  those  values  expected  for  Saturn's  dipole 
at  ~  5  Rs  to  10  Rs  [Smith  j^t  _al. ,  1980]  and  a  severe  decrease  of  ener¬ 
getic  proton  intensities  with  energies  >  100  keV  at  the  outer  edge  of  the 
plasma  torus  [Van  Allen  et  al. ,  1980a;  Simpson  et  al. ,  1980;  Tralnor  et 
al. ,  1980;  Filllus  et  al. ,  1980],  The  magnetic  field  depression  Is  pre¬ 
sumably  related  to  the  diamagnetic  effects  of  the  high-8  torus  plasma. 

For  example,  this  decrease  In  the  magnetic  field  intensity  is  ~  5  y  at 
8  R8  relative  to  the  expected  centered-dipole  total  field  of  41  y  [Smith 
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<5t  c»l. ,  1980].  Between  the  radial  distances  of  10  Rs  and  17  Rs  (magneto¬ 
pause)  the  measured  fields  exceed  the  dipole  values  by  ~  5  if.  A  detailed 
calculation  of  the  magnetic  perturbations  due  to  the  observed  ion  energy 
densities,  including  the  effects  of  the  magnetopause  current  system,  will  be 
undertaken  as  a  future  effort  in  order  to  show  that  the  measured  magnetic 
fields  are  consistent  in  detail  with  the  plasma  distributions  reported  here. 
Decreases  in  energetic  proton  intensities  are  not  as  easily  identified  with 
the  presence  of  the  plasma.  For  example,  the  loss  of  energetic  protons  may 
be  associated  with  strong  pitch-angle  scattering  into  Saturn's  ionosphere  by 
wave-particle  Interactions,  or  may  be  the  result  of  charge-exchange  with  an 
ambient,  diffuse  atomic  hydrogen  gas.  Observations  of  pitch-angle  distribu¬ 
tions  of  the  plasmas  and  of  plasma  waves  are  necessary  to  resolve  this  mech¬ 
anism  unambiguously.  Such  measurements  are  not  available  with  the  Pioneer 
11  spacecraft.  Wave-particle  interactions  should  be  abundant  in  such  high-3 
plasmas  if  the  earth's  plasma  sheet  is  used  as  a  guide  [cf.  Gurnett  et  al. , 
1976].  The  investigators  for  the  energetic  particle  instruments  cited  above 
also  report  remarkable  asymmetries  in  the  intensities  of  both  energetic  pro¬ 
tons  and  electrons  for  the  inbound  and  outbound  segments  of  the  Pioneer  11 
trajectory  at  radial  distances  ~  4  to  10  Rs  where  significant  azimuthal 
drift  asymmetries  and  large  effects  from  solar  wind  perturbations  are  not 
expected.  We  believe  that  this  asymmetry  is  probably  the  signature  of  a 
large  temporal  change  in  this  part  of  Saturn's  magnetosphere,  which  is 
associated  with  a  substantial  outflow  of  torus  plasmas  between  the  periods 
of  the  inbound  and  outbound  trajectories  as  the  solar  wind  pressure  on  the 
magnetosphere  decreased  [Wolfe,  et  aJL. ,  1980]  and  internal  magnetic  stresses 


lessened 


IB 


Rigid  corotation  of  the  plasma  allows  identification  of  the  M/Q  of  the 
dominant  species  within  the  torus:  ~  8  at  4.1  to  7.1  Rs,  and  ~  5  at  7.1  to 
8.2  Rs,  To  specifically  identify  the  positive  ions  we  are  guided  by  obser¬ 
vations  of  water  frost  or  ice  on  the  surfaces  of  Dione  and  Tethys  and  on 
the  ring  material  and  conclude  that  the  species  with  M/Q  ~  8  is  0^+,  and 
that  with  M/Q  ~  5  is  0^+.  Except  at  the  outer  edge  of  the  torus  at  ~  8.2 
to  8.6  Rs,  the  anticipated  ion  H*  from  photodissociation  of  H2O  was  not  ob¬ 
served  in  the  torus.  A  corresponding  upper  bound  on  the  H-*"  densities  in 
the  center  of  the  torus  is  ~  5  cm-^.  If  we  assume  that  the  relatively  large 
ion  thermal  speeds  corresponding  to  ion  temperatures  as  high  as  5  X  10^  °K 
do  not  severely  violate  the  conditions  for  Jordon's  [1969]  calculations  of 
ionization  equilibrium  for  oxygen,  the  dominant  charge  state  can  be  used  to 
obtain  a  gross  estimate  of  electron  temperature  at  the  orbits  of  Tethys  and 

Dione,  ~  4-8  X  10^  °K.  Similarly  the  dominance  of  0^+  at  ~  7  Rs  yields  an 

electron  temperature  of  ~  10^  °K.  These  electron  temperatures  are  also 
consistent  with  an  upper  bound  on  the  magnitude  of  the  spacecraft  potential 
relative  to  the  ambient  medium  of  ~  10  to  20  volts  as  indicated  by  the  good 
fits  of  the  observed  ion  velocity  distributions  at  energies  100  eV  to 
unperturbed  isotropic  Maxwellian  distributions. 

Our  current  interpretation  of  the  presence  of  a  large  oxygen  torus 
which  engulfs  the  orbits  of  Dione  and  Tethys  is  that  water  is  photodisso- 
clated  or  dissociated  by  charged  particles  as  a  source  of  atomic  hydrogen 

H  and  the  radical  OH.  These  products  then  are  ionized  by  either  electron 

impact  and/or  photoionization.  Overall  the  qualitative  nature  of  the 
temperature  and  density  profiles  in  Saturn's  oxygen  torus  are  similar  to 


those  for  the  sulfur  Ion  torus  within  Jupiter's  magnetosphere  at  Io's 
orbit  [cf.  Bagenal  et^  a_l . ,  1980].  In  particular  the  steep  temperature 
gradient  extending  from  approximately  Dione's  orbit  to  4.1  Rs  is  similar 
to  that  observed  in  the  vicinity  of  Io.  In  the  present  case  too,  the 
temperature  ratio  of  a  factor  of  10  for  temperatures  at  7  Rs  relative  to 
those  at  4  Rg  ts  much  greater  than  that  expected  for  simple  corotatlonal 
pickup  of  an  ionized  neutral.  The  corotatlonal  energies  for  oxygen  ions 
at  4  Rs  and  7.5  Rg  are  135  eV  and  470  eV,  respectively.  Bagenal  eit  al. 
[1980]  interpret  the  temperature  gradient  inside  the  orbit  of  Io  as  the 
inward  diffusion  and  radiative  cooling  of  sulfur  ions.  For  Saturn's  plas¬ 
ma  torus,  a  similar  interpretation  invokes  the  diffusion  and  radiative 
cooling  of  oxygen  ions  inward  from  the  orbit  of  Dione.  We  can  roughly 
estimate  the  order  of  magnitude  of  the  diffusive  loss  of  plasma  via  the 
relationship  for  the  diffusion  flux  F  ■  -D  AN/AR  where  D  is  the  diffusion 
coefficient,  and  AN/AR  Is  the  density  gradient.  We  take  AN  ~  50  cm~3  and 

AR  ~  3  Rs  from  our  present  observations  and  assume  that  D  ~  10~9  to 
2 

10  10  Rs/ sec.  This  latter  assumption  for  the  value  of  D  is  based  upon 
2 

values  ~  10"10  Rs/sec  for  protons  with  energies  >  80  MeV  in  Saturn's 
magnetosphere  [Van  Allen  et^  al. ,  1980b] ,  and  is  used  here  in  lieu  of  a 
more  satisfactory  means  of  estimating  D  [cf.  Siscoe,  1978],  The  order  of 
magnitude  of  the  diffusion  flux  F  becomes  3  X  10^  to  3  X  10^  ions/cm^-sec. 
If  we  assume  that  the  thickness  of  the  torus  is  4  Rs  then  the  total  ion 
flux  inwards  from  the  torus  toward  Saturn  is  ~  1.5  X  10^3  to  1.5  X  10^ 
ions/sec.  From  this  crude  estimate  for  diffusive  loss  we  can  estimate 
the  source  strength  at  the  surface  of  Dione.  Carlson  [1980]  has  consider¬ 
ed  various  mechanisms  for  the  dissociation  of  water  on  ice  surfaces  of 
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the  ring  material,  Including  photo-sputtering  [Harrison  and  Schoen,  1967], 
magnetospherlc  Ion  sputtering  [Cheng  and  Lanzerottl,  1978]  and  sublima¬ 
tion  with  subsequent  photodlssoclatlon  [Blamont,  1974],  and  concluded 
that  photo-sputtering  provides  the  largest  source  rate.  The  corresponding 
flux  escaping  from  the  water  frost  is  ~  10^  H  atoras/cm^-sec  for  a  molecu¬ 
lar  photodissociation  rate  of  1.1  X  10“7/molecule-sec.  If  we  assume  that 
(1)  a  similar  mechanism  is  operative  on  Dione,  (2)  the  OH  flux  is  equal 
to  that  for  atomic  hydrogen  and  (3)  the  diffusive  loss  rate  estimated 
above  is  the  dominant  loss  rate  for  the  plasma  torus,  then  the  flux  of 
OH  radicals  at  Dione' s  surface  must  be  ~  2  X  10?  to  2  X  10^/cm^-sec 
and  consistent  with  the  capabilities  of  the  photo-sputtering.  Since 
Tethys  is  characterized  with  a  similar  mass  and  radius  as  those  for  Dione 
[cf.  Cruikshank,  1978]  and  is  also  believed  to  be  ice-covered,  this  sat¬ 
ellite  would  be  expected  to  provide  a  source  of  oxygen  ions  of  similar 
magnitude.  It  is  noted  here  that  if  our  estimates  of  electron  thermal 
energies,  "  5  to  10  eV,  as  inferred  from  the  oxygen  charge  states  are 
correct,  then  the  surfaces  of  both  satellites  Dione  and  Tethys  are  exposed 
to  an  electron  flux  ~  10^/cm2_sec  with  these  energies  and  could  produce 
a  significant  contribution  to  the  molecular  dissociation  rates  via  electron 
impact.  The  binding  energy  for  H-OH  is  5.16  eV  [Dennefeld,  1974].  The 
radial  extent  of  the  plasma  torus  is  consistent  with  the  photodissociation 
of  water.  If,  for  example,  a  mechanical  energy  of  3  eV  is  available  after 
photodissociation  and  excitation  of  the  vibrational  and  electronic  states 
of  OH,  these  radicals  can  be  injected  into  elliptical  orbits  spanning  the 
radial  distance  range  of  ~  5.5  to  7.3  Rs.  Atomic  hydrogen  acquires  escape 
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velocity,  a  fact  that  can  account  for  the  general  relative  absence  of 
protons  in  the  plasma  torus.  On  the  other  hand,  In  the  simplest  model,  the 
Ion  temperature  due  to  corotatlonal  pickup  of  0+  upon  Ionization  of  the  OH 
radical  Is  expected  to  be  proportional  to  R^  as  contrasted  to  the  observed 
R6  variation  for  R  <  7  Rs. 

The  principal  reasons  for  considering  Dione  and  Tethys  as  primary 
sources  for  the  oxygen  torus  surrounding  Saturn  are  (1)  the  very 
suggestive  peaks  In  Ion  densities  at  their  orbits  (cf.  Figure  8),  (2) 
their  surfaces  are  known  to  be  covered  with  water  ice  or  frost  and  thus 
an  apparent  source  of  oxygen  via  dissociation  and  Ionization  and  (3)  the 
ion  density  and  temperature  radial  profiles  are  qualitatively  similar  to 
those  for  sulfur  ions  at  the  orbit  of  Io  in  Jupiter's  magnetosphere.  On 
the  other  hand,  the  temperature  profile  in  the  vicinity  of  the  orbits  of 
these  two  satellites  shows  no  clear  signature  of  their  presence.  Since 
Rhea  is  also  covered  with  water  frost  [cf.  Cruikshank,  1978]  and  its  mass 
and  radius  are  only  factors  of  2.1  and  1.6,  respectively,  greater  than 
those  of  Dione  [Anderson,  1980],  substantial  ion  densities  are  also 
expected  if  the  principal  source  Is  photo-sputtering  of  water  on  the 
satellite  surfaces.  No  such  signature  in  the  ion  densities  is  found  at 
Rhea's  orbit.  Sublimation,  a  source  of  lesser  magnitude,  Is  appreciably 
affected  by  the  radius  and  mass  of  the  body.  It  is  quite  possible,  in 
view  of  the  findings  of  substantial  atomic  hydrogen  densities  in  the 
vicinity  of  the  rings  [Judge  _et  al_. ,  1980],  that  Dione  and  Tethys  are 
only  minor  contributors  to  the  overall  oxygen  population  of  the  plasma 
torus.  Earth-based  measurements  of  the  resonantly-scattered  Lyman-a 
emissions  of  this  hydrogen  cloud  are  reported  by  Welser  et^  al.  [1977]. 


22 


In  this  Interpretation  the  rapid  density  gradient  at  4  Rg  in  the  vicinity 
of  the  orbit  of  Enceladus  is  due  to  a  combination  of  decreasing  ion 
temperatures  and  corotation  speeds  to  the  extent  that  the  ion  ener¬ 
gies  are  below  the  energy  threshold,  E/Q  =  100  eV,  of  the  plasma  Instru¬ 
ment.  This  possibility  is  strikingly  suggested  in  the  character  of  the 
spectrograms  of  Plate  1  for  ~  1400  ERT.  The  outer  edge  of  the  A  ring  is 
positioned  only  about  1.5  Rs  away  from  the  termination  of  usable  respons¬ 
es  of  the  plasma  instrument  at  4.1  Rg.  The  source  rate  necessary  to  sup¬ 
port  the  atomic  hydrogen  densities  in  the  vicinity  of  the  rings  is  ~  3  X 
10^8  atoms/sec,  and  is  consistent  with  the  mechanism  of  photo-sputtering 
of  ring  ice  material  [Carlson,  1980],  This  source  rate  is  a  factor  of  ~ 

3  X  104  greater  than  that  from  photo-sputtering  from  the  surface  of  Dione 
or  Tethys.  The  actual  number  of  oxygen  ions  from  ionization  of  OH  and 
escaping  impact  with  ring  material  is  difficult  to  estimate  without 
further  calculations.  In  this  interpretation  the  ions  diffuse  outward 
from  the  rings,  are  rapidly  heated  as  a  function  of  increasing  radial 
distance  and  subsequently  gain  sufficient  thermal  and  corotational 
energies  to  be  observed  with  the  plasma  instrument  at  ~  4  Rs.  Thus 
heating  of  the  oxygen  ions  is  required  as  the  ions  diffuse  radially 
outward  to  ~  6.5  Rg  where  a  high-8  plasma  occurs. 

In  summary,  due  to  (1)  the  monotonically  increasing  temperatures 
with  increasing  radial  distance  in  the  torus,  which  are  unaffected  by  the 
presence  of  Dione  and  Tethys  and  (2)  the  substantially  greater  potential 
source  rate  for  OH  radicals  at  the  rings,  we  favor  the  rings  as  the  pri¬ 
mary  source  of  torus  plasmas  and  suggest  Dione  and  Tethys  as  secondary 


sources  as  evidenced  by  increases  of  oxygen  ion  densities  at  their  or¬ 
bits.  This  interpretation  will  also  yield  the  substantially  lesser  oxy¬ 
gen  ion  temperatures  that  are  observed  relative  to  those  expected  for 


corotational  pickup  since,  at  the  outer  edge  of  the  A  ring  at  2.3  Rs,  the 
neutral  particle  orbital  speed  differs  from  the  corotation  speed  by  only 
~  7  km/sec  which  corresponds  to  an  oxygen  ion  temperature  of  ~  5  X  10^ 

°K.  These  observations  indicate  the  possibility  of  a  torus  of  H2O  and  OH 
in  the  vicinity  of  the  orbits  of  Dione  and  Tethys  with  substantially  les¬ 
ser  densities  than  those  at  the  rings.  The  diffuse  hot  component  of  the 
oxygen  velocity  distributions  detected  with  the  plasma  instrument  in  the 
torus  (cf.  Figure  7)  could  be  the  signature  of  the  ionization  of  OH  from 
Dione  and  Tethys.  For  example  the  oxygen  temperature  corresponding  to 
corotational  pickup  with  ionization  of  OH  is  ~  3  X  10^  °K.  The  observed 
temperature,  if  the  ions  in  the  high-speed  tail  of  the  velocity  distribu¬ 
tion  are  o2+,  is  ~  5  X  106  °K. 

The  source  of  the  protons  observed  at  radial  distances  beyond  8.8  Rs 
is  not  clearly  resolved  here.  One  possibility  is  that  these  protons  are 
the  ionization  products  of  a  tenuous  gas  of  atomic  hydrogen  from  the  dis¬ 
sociation  of  water  frost  on  the  ring  material.  Photo-sputtering  will 
allow  atomic  hydrogen  to  gain  escape  velocity  and  thus  provide  a  source 
mechanism  at  these  radial  distances.  The  substantially  fluctuating  ion 
temperatures  and  densities  at  R  >  8.8  Rg  suggest  rapid  temporal  varia¬ 
tions  in  the  topology  of  the  outer  magnetosphere.  Other  possible  sources 
for  these  low-energy  protons  include  the  magnetosheath,  Titan  (the  day- 
side  magnetopause  was  crossed  inside  Titan's  orbit  during  a  period  of  en¬ 
hanced  solar  wind  pressure),  and  Saturn's  ionosphere.  With  the  current 
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measurements,  it  does  not  appear  possible  to  clearly  distinguish  among 
these  four  sources  for  ions  in  the  outer  magnetosphere. 

In  the  spirit  of  speculation,  it  is  noted  here  that  if  the  plasma 
torus  extends  inward  to  the  rings  of  Saturn  the  satellite  Enceladus  will 
be  nonuniformly  exposed  to  a  large  flux  of  oxygen  ions  on  its  trailing 
face  since  the  corotational  bulk  flow  speed  exceeds  the  satellite  orbital 
speed  by  28  km/sec.  There  is  evidence  that  the  trailing  face  of  Encela- 
dus  is  brighter  due  to  the  presence  of  water  frost  than  the  leading  face 
[Cruikshank,  1978].  The  oxygen  flux  is  estimated  here  at  ~  1.5  X  10^ 
ions/cm^-sec.  For  comparison  the  photo-sputtering  loss  rates  are  ~  1  X 
10**  molecules/cm^-sec  and  the  sublimation  rate  for  water  frost  at  100°K 
is  ~  5  X  10^  molecules/cm^-sec  [Dennefeld,  1974],  It  is  possible  that 
the  recombination  of  the  oxygen  and  hydrogen  produces  water  frost 
preferentially  on  the  trailing  face. 

A  constant  source  of  oxygen  ions  and  their  subsequent  energization 
in  the  plasma  torus,  and  the  observational  fact  that  the  8  of  the  plasma 
torus  is  of  the  order  of  unity  as  deep  in  Saturn's  magnetosphere  as  6.5 
Rs,  suggest  that  the  plasma  pressure  may  be  periodically  relieved  by 
spontaneous  plasma  outflow  through  the  magnetotail.  We  note  here  that 
the  condition  for  plasma  outflow  in  a  rotating  dipole  field  [Hill  et  al. , 
1974]  is  approximately  met  over  the  radial  distance  range  6  to  8  Rs. 

This  condition  requires  that  the  corotational  energy  density  of  the  plas¬ 
ma  is  equal  to  or  greater  than  the  energy  density  of  the  dipole  magnetic 
field.  Major  reconfigurations  of  the  magnetosphere  beyond  these  radial 
distances  can  also  be  expected.  The  total  number  of  oxygen  ions  in  the 


torus,  assuming  a  thickness  of  4  Rs  is  ~  2  X  1033  ions.  If  the  source  of 

ions  is  photo-sputtering  from  the  surface  of  Dione  and  in  the  unlikely 

situation  that  there  are  no  significant  losses,  the  torus-filling  time  is 

~  (2  X  1033  ions)/(7  X  10^3  ions/sec)  5  3  X  10^  seconds  (100  years).  On 

the  other  hand,  if  the  oxygen  source  rate  is  only  17,  of  the  source  rate 

for  atomic  hydrogen  at  the  rings  as  estimated  by  Carlson  [1980],  then 

this  torus-filling  time  is  ~  (2  X  1033  ions)/(3  X  10^6  ions/sec)  *  7  X 

10^  seconds  (81  days).  The  corresponding  diffusion  coefficient  D  is  ~  5 
7  2 

X  10“'  Rs/sec  and  the  average  radial  drift  speed  is  ~  10“*  km/sec.  If 
the  absence  of  density  minima  at  the  satellite  orbits  is  interpreted  in 
terras  of  drift  speeds  sufficiently  great  to  provide  displacements  com¬ 
parable  to  or  greater  than  the  diameter  of  Dione  or  Tethys  in  one  plasma 
corotational  period,  the  corresponding  lower  limit  on  the  drift  speed  Is 
~  (1000  km)/(3.7  X  10^  sec)  *  3  X  10“  ^  km/sec.  Thus  a  ring  source  for 
the  oxygen  ions  can  profoundly  affect  raagnetospheric  configuration  and 
dynamics  beyond  ~  6  Rs.  The  next  series  of  situ  measurements  will  be 
gained  with  the  Voyager  spacecraft.  The  above  calculations  indicate  that 
the  torus  configuration  could  be  greatly  different  relative  to  that  re¬ 
ported  above  for  the  Pioneer  11  encounter. 


Appendix:  Plasma  Analyzer  Responses  to  a 
Corotating  Maxwellian  Velocity  Distribution 

We  summarize  here  our  treatment  of  the  responses  of  the  hlgh- 
resolution  electrostatic  analyzer  to  a  corotating  plasma  within  Saturn's 
magnetosphere.  For  the  measurements  presented  herein  the  center  of  the 
f ield-of-view  of  sensor  1  Is  directed  at  50.4°  to  the  spacecraft  spin 
axis.  The  responses  of  this  sensor  are  accumulated  onboard  only  for  each 
of  two  one-half  spin  periods  (see  Figure  1).  The  average  directional, 
differential  Intensity  of  positive  Ions  for  a  half  spin  period,  dJ/dEdfi, 
is  related  to  the  corresponding  accumulated  counts,  C, 

dj  «  4.64X10*C  ,  (1) 

dEdSJ  w2 

where  M  and  V  are  the  ion  mass  and  velocity,  respectively.  Cgs  units  are 
employed  throughout  the  discussion.  The  average  phase  space  density, 
n(V)  In  units  of  sec^/cm^,  is  related  to  above  intensities. 


n(V)  "  M  dJ  » 
V2  dEdft 

thus  yielding 


(2) 


n(V) 


4.64X10*C 

V4 


O) 
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For  a  given  V,  or  E/Q,  two  determinations  of  the  average  n(V)  are 
observed  with  the  plasma  analyzer,  one  each  for  the  two  half-spin 
accumulation  periods. 

We  must  now  find  the  expected  responses  of  the  instrument  to  a 
corotating  Maxwellian  plasma  with  number  density  N,  temperature  T  and 
corotating  velocity  vector  in  the  spacecraft  reference  frame.  This 
corotation  speed  is  sufficiently  large  to  allow  neglect  of  the  spacecraft 
velocity  vector  relative  to  Saturn  in  the  present  approximation.  For  the 
corotating  Maxwellian 

n(V)  “  N  exp(-MV,2/2kT)  ,  (4) 

^2skT  j3/2 

where  *  ^-V0.  Note  that  for  a  given  value  of  n(V),  the  temperature  T 
is  proportional  to  mass  M.  Hence  two  ions  with  the  same  M/Q  value,  for 
example  H2+  and  He2+,  will  have  identical  velocity  distributions  if 
T(He2+)  «*  2T(H2+).  Thus  the  ion  must  be  identified  for  an  accurate 
assessment  of  the  temperature  T.  Since  the  electrostatic  analyzer  is 
capable  of  assessing  only  M/Q,  determination  of  the  plasma  temperature  T 
depends  critically  on  the  mass  of  the  assumed  ion.  The  average  densi¬ 
ties  n(V)  for  the  two  half-spins  are 


n(V)  - 


'+1T 

n(V)d$ 


(5) 
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with  »  n/2  and  3x/2,  respectively.  The  angle  $  *8  the  6pin  phase 
angle  for  the  spacecraft  rotation  and  is  referenced  to  <t>  ■0®  when  the 
sensor  f ield-of-view  is  parallel  to  the  ecliptic  and  ascending  toward  the 
north  ecliptic  pole. 

In  order  to  numerically  integrate  equation  (5),  ^o*V  must  be  determined 
as  a  function  of  spin  phase  angle  The  ecliptic-equinox  coordinate  system 
(epoch  1950.0)  can  be  conveniently  utilized  to  evaluate  The  Xe-axis 

is  directed  toward  the  vernal  equinox  and  the  Ze-axis  toward  the  north 
ecliptic  pole.  The  positive-ion  velocity  vectors  Ve  viewed  by  the  sensor  as 
a  function  of  $  in  these  coordinates  are 

-V  cos8q  cosB  cosa  -  V  sin0o  cos^  sinB  +  V  sin0o  sin<J>  cosB  sina 
Ve  «  V  cos0o  sinB  cosa  -  V  sin0o  cos$  cosB  -  V  sin0o  simj>  sinB  sina  (6) 

-V  cos0o  sinaQ  -  V  sin0o  sin$  cosa 

with  0O  «*  50.4®  and  spacecraft  spin  axis  direction  a  «  -1.8®,  B  *  14.7®.  In 

order  to  find  the  corotational  velocity  vector  VQ  in  the  ecliptic-equinox 

coordinates  it  is  necessary  to  first  find  the  spacecraft  position  in  Saturn 

spin-axis  coordinates,  Xs  ,  then  to  determine  in  these  coordinates  and 

sc  o 

subsequently  transform  this  vector  into  ecliptic-equinox  coordinates,  ^e. 

The  spacecraft  position  is 

cosOe  cos<J>e  Xe  +  cosOe  sin$e  Ye  -  sin0e  Ze 
sc  sc  sc 

^s  m  -8in<j>e  Xe  +  cos$e  Ye 

SC  SC  8C 

sinOe  cos$G  Xe  +  sin0e  sln<|>e  Ye  +  cos0e  Ze 
sc  SC  8C 


«/*■  ia.  mV.* 


<7) 
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where  Qe  =  28.1°,  $e  •  78.8°  is  the  direction  of  the  Saturn  spin  axis  and 

is  the  spacecraft  position  in  Saturn-centered  ecliptic-equinox 
sc 

coordinates.  The  corotational  velocity  in  ecliptic-equinox  coordinates 
becomes 


-wR 

cos6  8 

siw}>s 

COSlj>SS 

cos6ss  + 

coR 

cosOs 

cosij)8 

sin<j>ss 

sc 

sc 

sc 

sc 

?e  - 

wR 

cos6s 

sin<j>s 

sin<j)SS 

cosO88  + 

ioR 

o 

o 

03 

CD 

03 

COSlfl3 

COS(J>SS 

o 

sc 

sc 

sc 

sc 

-o>R 

cos68 

sin<|>s 

sin9ss 

sc 

sc 

(8) 


where  the  Saturn  angular  rotation  speed  w  «  1.73  X  10“4  rad/sec,  R  is  the 

radial  distance  from  Saturn  to  the  spacecraft,  9s  =  Arcsin  (Zs  /R),  <|>3  « 

sc  sc  sc 

Arctan  (Ys  /Xs)  and  the  Saturn  spin  axis  direction  9SS  =*  -28.1°, 
sc  sc 

*ss  «*  -78.8°. 

The  above  transformations  (6)  and  (8)  allow  numerical  integration  of 
the  two  integrals  (5)  which  are  directly  related  to  the  observed  average 
densities  (3).  Increments  of  A<j>  =  1.2  X  10“2  radian  are  chosen  for  the 
numerical  integration.  A  nonlinear  least  square  fit  which  utilizes  the 
method  proposed  by  Gauss  (courtesy  of  K.  Sando,  Chemistry  Department)  [cf. 
McWilliams  et  al.,  1965]  is  employed  to  determine  N  and  T  from  the  integral 
(5)  corresponding  to  the  angular  sector  closest  to  the  direction  of 
corotation,  <j>'  »  x/2,  for  a  given  ion  mass  M.  This  best-fit  computation 
for  N  and  T  is  repeated  for  M  ■  1,  2,  4,  8  and  16  amu.  The  plasma  analyzer 
is  not  capable  of  direct  identification  of  the  ion  species.  However,  the 
Information  concerning  angular  distributions  available  from  the  two  inte¬ 
grals  (5),  <j>’  -  it/2  and  ■  3n/2,  with  the  above  assumption  of  a  coro¬ 
tating  Maxwellian  distribution,  allows  an  indirect  determination  of  the 


mass-per-unit  charge  (M/Q)  values  for  the  positive  Ions.  Hence  the  ob¬ 
served  n(V),  (3)  above,  In  the  spin  sector  generally  directed  away  from  the 
corotation  can  be  used  for  comparison  with  the  computed  n(V)  for  4>'  “  3ir/2 
to  find  an  approximate  value  for  M/Q  by  imposing  the  best-fit  N  and  T  from 
(5)  for  <J>'  =  ir/2.  The  counting  statistics  for  measurements  corresponding 
to  the  integral  <J>'  *  3it/2  are  not  always  sufficient  for  such  species 
identification  (cf.  Figure  6).  These  extensive  numerical  computations  are 
performed  for  each  of  the  62  individual  measurements  of  ion  velocity  dis¬ 
tributions  within  Saturn’s  dayside  magnetosphere  as  presented  In  the  text. 
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Plate  Caption 

Plate  1.  Color-coded  energy-time  spectrograms  for  the  responses  of  the 
hlgh-resolutlon  plasma  analyzer  on  board  Pioneer  11  to  low- 
energy  positive  Ions  as  this  spacecraft  traverses  the  dayside 
magnetosphere  of  Saturn  during  0000  to  1650  ERT  (0T  of  re¬ 
ceipt  of  signal  at  earth)  on  September  1,  1979.  (See  text 
for  details  of  spectrograms.)  The  upper  spectrogram  sum¬ 
marizes  responses  during  which  the  f ield-of-view  for  sensor  1 
of  the  plasma  analyzer  is  generally  directed  toward  the  coro- 
tational  flow,  and  the  bottom  spectrogram  displays  responses 
when  the  f ield-of-view  is  scanning  directions  at  larger  an¬ 
gles  to  the  corotational  flow  vector.  A  large  torus  of  oxy¬ 
gen  ions  is  encountered  during  the  period  ~  1000  to  1500  ERT. 
The  high  responses  at  ~  1500  to  1630  ERT  are  due  to  penetrat¬ 
ing  energetic  charged  particles.  The  responses  in  the  upper 
spectrogram  are  identified  in  the  present  analysis  as  those 
corresponding  to  angular  sector  =  ir/2,  and  those  of  the 
bottom  spectrogram  to  angular  sector  <J>'  *  3n/2. 


Note:  This  plate  is  to  be  published  in  color. 
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The  left-hand  panel  shows  the  trajectory  of  Pioneer  11  as 
projected  onto  Saturn's  equatorial  plane  during  the  period  of 
plasma  measurements  displayed  In  Plate  1.  Angular  sampling 
of  the  plasma  velocity  distributions  with  sensor  1  of  the 
plasma  analyzer  Is  summarized  with  the  right-hand  diagrams. 
The  axis  of  the  sensor's  f leld-of-view  is  directed  at  an 
angle  50.4°  to  the  spacecraft  spin  axis.  The  responses  of 
the  sensor  are  accumulated  for  two  angular  sectors  each 
corresponding  to  one-half  spin  period  of  the  spacecraft  and 
referenced  to  the  ascending  crossing  of  the  ecliptic  plane. 
These  two  angular  sectors  are  referred  to  herein  as  <J>’  =  x/2 
and  *  3*/2,  respectively. 

Best  nonlinear  least-squares  fits  of  corotating  hydrogen 
velocity  distributions  with  various  fractions  of  the  rigid 
corotation  speed  VQ  to  the  observed  Ion  velocity  distribu¬ 
tions  in  the  outer  magnetosphere  at  15.7  Rs  (1  Rs  *  Saturn 
radius  =  6  X  10^  km).  The  corotation  speed  V0  Is  160  km/sec. 
A  good  fit  Is  obtained  for  a  speed  of  0.5  VD>  or  80  km/sec. 
Continuation  of  Figure  2  for  the  ion  velocity  distribution 
observed  at  0117  ERT  (15.1  Rs).  A  good  fit  to  a  corotating 
Maxwellian  is  not  found  for  this  series  of  observations. 
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Best  nonlinear  least-squares  fits  for  the  ion  velocity 
distributions  observed  at  7.1  Rs  in  the  plasma  torus, 
assuming  values  of  M/Q  =  1(H+),  4(He+),  8(f)2+)  and 
16(0+).  The  overall  best  fit  is  for  a  species  with  M/Q 
between  4  and  8. 

Continuation  of  Figure  4  for  ion  velocity  distributions 
observed  at  1219  ERT  (6.4  Rs).  The  best  overall  fit  is 
for  ()2+  with  a  density  N  =  40.8  cm-^  and  a  temperature  T 
*  1.6  X  106  °k.  The  corotation  speed  is  65  km/sec. 
Continuation  of  Figure  4  for  five  representative  ion 
species  at  5.6  Rs  where  the  ion  temperature  is  suffi¬ 
ciently  low  that  no  usable  responses  are  available  f'  • 
angular  sector  <f> ’  ■  3n/2  which  views  large  angles 
relative  to  the  corotation  bulk  flow. 

Best  fits  for  the  high-speed  tails  of  the  ion  velocity 
distribution  at  6.6  Rs.  These  ions  are  most  likely  0^+ 
or  0^+  corresponding  to  the  range  of  best-fit  M  values 
of  ~  4  to  8  amu.  A  substantial  contribution  to  this  ion 
velocity  distribution  by  H+  or  0+  ions  appears  to  be 
excluded. 


Ion  densities,  temperatures  and  dominant  species  as 
functions  of  planetocentric  radial  distance.  Identifi¬ 
cation  of  the  dominant  ion  species  is  discussed  in  de¬ 
tail  in  the  text.  The  oxygen  torus  extends  from  ~  4  Rg 
to  8  Rs.  The  rapid  decline  of  densities  at  ~  4  Rs  may 
be  due  to  ion  energies  decreasing  with  decreasing  radial 
distance  to  values  below  the  energy  (E/Q)  threshold  of 
the  plasma  instrument  at  100  eV-per-unit  charge.  Ion 
temperatures  within  the  plasma  torus  at  ~  4  Rg  to  7  Rs 
monotonically  decrease  with  decreasing  radial  distance. 
The  outer  edge  of  the  A  ring  is  located  at  2.3  Rs. 

The  ion  energy  density,  e  in  units  of  ergs/cm^,  and  the 
ratio  of  plasma  and  magnetic  field  pressures,  0,  as 
functions  of  radial  distance.  Note  the  relatively  high 
values  of  0  in  the  plasma  torus  as  close  to  the  planet 
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